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The lattice cluster theory (LCT) for the thermodynamics of polymer systems has recently been reformulated 
to treat strongly interacting self-assembling polymers composed of fully flexible linear telechelic chains [J. 
Dudowicz and K. F. Freed, J. Chem. Phys. 136 , 064902 (2012)]. Here, we further extend the LCT for linear 
telechelic polymer melts to include a description of chain semiflexibility, which is treated by introducing 
a bending energy penalty whenever a pair of consecutive bonds from a single chain lies along orthogonal 
directions. An analytical expression for the Helmholtz free energy is derived for the model of semiflexible 
linear telechelic polymer melts. The extension provides a theoretical tool for investigating the influence of 
chain stiffness on the thermodynamics of self-assembling telechelic polymers, and for further exploring the 
influence of self-assembly on glass formation in such systems. 


I. INTRODUCTION 

Telechelic polymers, containing one bifunctional asso¬ 
ciative group at each of the chain ends (called “stick¬ 
ers” ), provide a striking example of macromolecules that 
are capable of supramolecular self-assemblyji The self- 
assembly in telechelics has attracted considerable atten¬ 
tion for both technological and fundamental reasons. On 
one hand, the distinctive characteristics of telechelics, 
arising from the reversible formation and breakage bonds 
during the dynamical self-assembly, open the prospect of 
many new applications)^^— which are generally not acces¬ 
sible by conventional polymerization. On the other hand, 
the theoretical description of the self-assembly process in 
telechelics represents a large challenge because of the in¬ 
terplay of the strong interactions between the stickers 
and the van der Waals interactions between the polymer 
segments. Additional complexity arises from the internal 
chemical structure of individual telechelic molecules as 
well as from the additional reversible chain connectivity 
introduced by the associative clusters. 

Theories^^— and simulationsi^^— of self-assembly in 
telechelic polymers traditionally employ highly coarse 
grained models that represent the assembling molecular 
species as a structureless entity and hence, neglect the in¬ 
fluence of local molecular structure on the self-assembly 
of telechelic polymers. However, a deep understanding 
of the relation between molecular structure and physi¬ 
cal properties is important in guiding the rational design 
of telechelic polymer materials. To address this need, 
Dudowicz and Freed^l have recently developed an inter¬ 
mediate level of coarse grained models that retain min¬ 
imal aspects of molecular structure and interactions in 
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telechelic polymers. The initial implementation of the 
model has proceeded by extending the lattice cluster 
theory (LCT) for the thermodynamics of polymer sys- 
tems^Sr.^ to strongly interacting, self-assembling poly¬ 
mers composed of fully flexible linear telechelic chains. 
The description provided by this extended LCT enables 
establishing the relation between the molecular structure 
dependent interaction parameters of the model and the 
thermodynamic properties of these complex fluidsi ^^i^^ 
Therefore, the extended LCT for telechelic polymers fills 
a gap between the limited predictive abilities of the pre¬ 
vious theories with structureless monomers and the ther¬ 
modynamic complexity of the molecular details of the 
self-assembly processes. 

Several improvements are desirable within the LCT for 
telechelic polymers since the initial study by Dudowicz 
and Freed^i only considers the simplest case, i.e., fully 
flexible linear telechelic polymer chains. For example, the 
applicability of the theory will be enhanced by including 
the description of monomer structures, a description that 
has been demonstrated to greatly affect the thermody¬ 
namics of polymer systems— i Another improvement rep¬ 
resents the focus of the present paper, namely, allowing 
the polymer chains to be semiflexible. Our motivation, 
in part, comes from the fact that chain semiflexibility 
plays an important role in determining the self-assembly 
process and structural properties of telechelic polymers; 
e.g., computer simulationsii2^i2& indicate that the struc¬ 
ture of the aggregates in telechelic polymers changes sig¬ 
nificantly with the degree of chain stiffness. More im¬ 
portantly, chain semiflexibility is crucial for exploring 
glass formation in systems of self-assembling telechelic 
polymers within the generalized entropy theory (GET))2i 
which is a merger of the LCT for the thermodynamics of 
semiflexible polymers^ and the Adam-Gibbs (AG) re¬ 
lation between the structural relaxation time and the 
configurational entropy ;^^!^^ because the characteristic 
glassy behavior does not appear in the GET when the 
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FIG. 1. Schematic illustration of the lattice model of a single 
linear telechelic polymer chain with M united atom groups. 
Solid circles designate the chain’s ends that can participate in 
sticky interactions. Solid lines represent the bonds between 
united atom groups, while dotted line indicates the presence 
of intervening united atom groups and bonds between the two 
separated portions on the chain. The model prescribes two 
different nearest neighbor interaction energies e and Ss for 
the weakly interacting interior united atom groups and the 
stickers, respectively. 


polymer chains are modeled as being fully flexible, as 
noted in Ref. Therefore, the present paper further 
extends the LCT to describe the thermodynamics of a 
melt of semiflexible linear telechelic chains. The exten¬ 
sion in the present paper, thus, also enables investigat¬ 
ing the influence of self-assembly on glass formation in 
telechelic polymer melts in the GET. 


II. LATTICE MODEL OF SEMIFLEXIBLE LINEAR 
TELECHELIC POLYMER MELTS 

As in previous work,— the polymer chains are 
placed on a d-dimensional hypercubic lattice with Ni 
lattice sites, each with z = 2d nearest neighbors. The 
present work considers a compressible (i.e., the lattice 
contains empty sites) telechelic polymer melt, where the 
system consists of m semiflexible linear chains with M 
united atom groups per chain, producing the volume frac¬ 
tion of the polymer chains as (j) = mM/Ni. Semiflexibil¬ 
ity is treated by introducing a bending energy penalty 
Eh whenever a pair of consecutive bonds from a single 
chain lies along orthogonal directions.— Telechelic chains 
are modeled by distinguishing each chain’s end segments 
(called “stickers” and represented as solid circles in Fig. 
1) from those lying in the chain interior (depicted by 
open circles in Fig. 1), as introduced in Ref. More 
specifically, two stickers on nearest neighbor lattice sites 
can form a sticky “bond” with an attractive sticky inter¬ 
action energy that may greatly exceed the microscopic 
van der Waals interaction energy e, which describes the 
interaction strength between two non-stickers as well as 
between a sticker and a non-sticker (Fig. 1). The stick¬ 
ers are assumed to be mono-functional and hence, each 
sticker can only participate in one sticky interaction. In 
accord with the previous worl(— and the analysis of 
Jacobson and Stockmayer;^ both cyclic and linear as¬ 


sociative clusters may form. Because the lattice model 
is obviously coarse grained, the sticky energy is gener¬ 
ally a temperature dependent parameter, but taken as a 
constant for simplicity. An attractive microscopic sticky 
interaction energy Cg is defined here as negativej^l while 
the attractive microscopic van der Waals interaction en¬ 
ergy e is treated as positive, in line with the original 
LCT— Since the sticky bonds (which are essentially 
non-covalent) formed by stickers differ from the covalent 
bonds within a single chain in real telechelic polymers, 
semiflexibility constraints are not applied for a pair of se¬ 
quential bonds containing one sticky bond in the present 
work, in line with computer simulations of semiflexible 
telechelic polymers— When the sticky interactions 
are strong, the clusters formed by rigid rod telechelic 
chains become freely hinged chains that are connected by 
flexible joints. Therefore, crystallization may be avoided 
in the models of rigid linear telechelic polymers, a feature 
that is potentially important for exploring glass forma¬ 
tion in such systems. 

While the previous study of fully flexible linear 
telechelic chains with the LCT considers the system as 
incompressible solutions;^ the mathematical equivalence 
between a compressible melt and an incompressible so¬ 
lution allows drawing conclusions for both types of sys¬ 
tems. The model for an incompressible solution replaces 
the empty sites by one-bead solvent molecules instead. 
The excess free energy of a compressible polymer melt 
is isomorphic to that of an incompressible polymer solu¬ 
tion, with the e parameter being replaced by the exchange 
energy Cex = epp + ^ss-^Cps, where Cpp, egg and epg repre¬ 
sent the strengths of the nearest neighbor interaction be¬ 
tween two polymer segments, two solvent molecules and 
a polymer segment and a solvent molecule, respectively. 
Because the theory in the present paper will be used in 
a subsequent paper to explore glass formation in linear 
telechelic polymer melts, the discussion of the model and 
results are expressed in the present work with reference 
to compressible melts. 


III. LATTICE CLUSTER THEORY FOR SEMIFLEXIBLE 
LINEAR TELECHELIC POLYMER MELTS 

As demonstrated in a previous paper the develop¬ 
ment of the LCT for associating polymers begins by not¬ 
ing that the specific Helmholtz free energy / = F/Ni 
(with F the total Helmholtz free energy ) for a self¬ 
assembling system is the sum of the specific free energy 
fo in the absence of strong sticky interactions and a con¬ 
tribution fg arising from the sticky interactions, 

f = fo + fs. ( 1 ) 

Specializing to the model in the present work, fo is the 
specific free energy for a polymer melt of semiflexible 
linear chains, where each chain’s ends are indistinguish¬ 
able from the other united atom groups. This is our 
reference system. The free energy has been derived for 
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the model of multicomponent systems with structured 
monomer chains in the previous paper— with some cor¬ 
rections given in Ref. l38l . Therefore, this section first 
summarizes the result for the free energy fo specialized 
to the model of polymer melts composed of semiflexible 
linear chains, and then describes the evaluation of fs with 
the emphasis on how to include chain semiflexibility in 
the presence of sticky interactions. 


A. Free energy for the reference system 

The LCT yields the specific Helmholtz free energy fo 
of a semiflexible polymer melt in the general form,— >2^ 

6 

= ( 2 ) 

i=l 

where /3 = l/iksT) with ks being Boltzmann’s constant 
and T designating the absolute temperature. The term 
Pjmf represents the zeroth-order mean-field contribu¬ 
tion, and has the following form for a melt of semiflexible 
linear chains, 

(3) 

where N 2 is the number of runs of two consecutive bonds 
in a single chain, and Zb = (zp — 1) exp(—-|- 1 with 
Zp = zj^ and being the bending energy. The coeffi¬ 
cients Ci {i = 1,...,6) are obtained by collecting terms 
corresponding to a given power of <j), and these coef¬ 
ficients are generally a function of z, T, e, Eb, and a 
set of counting indices Ui = Ni/M {i = 1,...,4), where 
the counting factor Ni denotes the number of runs of 
i consecutive bonds in a single chain and simplifies to 
Ni = M — i for linear chains. The explicit expressions 
for Ci (i = 1,..., 6) for a melt of semiflexible linear chains 
are provided in Appendix A. 


B. LCT partition function for semiflexible linear telechelic 
polymer melts and its diagrammatic representation 

As explained in Refs. the main idea behind the 

extension of the LCT to strongly interacting polymer sys¬ 
tems lies in the extraction of terms associated with the 
strong interactions from the cluster expansion. Conse¬ 
quently, the partition function W for the self-assembling 
systems is first defined for systems with a constant num¬ 
ber iJ of sticky bonds and then is summed over all pos¬ 
sible El to yield W. 

Imposing the constraints associated with the pres¬ 
ence of excluded volume interactions, chain connectivity, 
semiflexibility, nearest neighbor van der Waals interac¬ 
tions, and sticky interactions, the LCT partition func¬ 
tion W{E[, Ns,m, M) for a melt of semiflexible linear 


telechelic chains is derived as 
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where the symbol X]{r-} represents a restricted sum over 
the positions of all united atom groups subject to the 
constraint prohibiting multiple occupancy of any lat¬ 
tice site, S is the Kronecker delta function (the physi¬ 
cal mea ning of each specific S is explained in detail in 
Refs, ill a is the lattice vector, the factor E is de¬ 
fined by A = exp(-/3L;b), /pair = exp(,5e) - 1 is the 
Mayer /-function. Ns is the total number of stickers in 
the system, is the sticker-sticker interaction energy, the 
index t labels the sticky bonds, and and are the 
positions of the nearest neighbor stickers p and g, which 
belong to the set {r^} of lattice sites that are occupied 
by the polymer chains. The factors 2“™ and 1/m! in Eq. 
(4) account for the indistinguishability of the two chain 
ends and of the chains themselves, while the combinato¬ 
rial factor Ns\/[{Ns — 2H)\E[\2^\ in Eq. (4) accounts for 
the number of ways of forming a total of H sticky bonds 
from the total number Ns of stickers in the system. 

In analogy to the treatment of the partition func¬ 
tion for the reference system^ the partition function 
W{H, Ns, m, M) for the telechelic polymers can schemat¬ 
ically be rewritten in the convenient form of multiple clus¬ 
ter expansions, 
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where the zeroth-order mean-field terms and 

are provided by Eq. (19) in Ref. HO and Eq. (11) in 
Ref. [ 23 , respectively. The terms A'bond, G, Afpair, K, 
kbend, and X^ticky in Eq. (11) are given by 
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FIG. 2. Contributing cumulant diagrams with a single sticky 
bond for the model of semiflexible linear telechelic polymer 
melts. Each diagram contains at least one bending constraint 
line but no interaction lines. Open and solid circles corre¬ 
spond to non-sticky and sticky united atom groups, respec¬ 
tively. Solid straight lines denote non-sticky or sticky bonds, 
while dotted straight lines in j and h indicate the presence 
of one or more intervening bonds within a chain between the 
depicted bonds. The diagrams are derived from Fig. 5 in 
Ref. [13. 
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FIG. 3. Coirtributing cumulant diagram with two sticky 
bonds, one bending constraint line, but no interaction lines 
for the model of semiflexible linear telechel ic p olymer melts. 
The diagram is derived from Fig. 6 in Ref. 1^ . 
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Equations (11-17) are achieved by taking advantage of 
the fact that each Kronecker delta function in Eq. (4) 
can be expressed in the convenient form, 

S = A+{6-A) = a(i+^~'^ 
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( 12 ) 


where the term A represents an average of the contribu¬ 
tion from S and thus defines a factor appearing in the 
zeroth-order mean-field term or 1E™1. 

Expanding the product in Eq. (5) leads natu¬ 
rally to the cluster expansion for the partition function 
WiH,N,,m,M), 
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FIG. 4. Contributing cumulant diagrams with a single sticky 
bond for the model of semiflexible linear telechelic polymer 
melts. Each diagram contains one bending constraint line 
and one interaction line. The diagrams are derived from Fig. 
8 in Ref. [H. 


where {o. — J, ...) is either Xuond; f^^pairj Uliaend; 
or Xgticky The systematic calculation of contributions 
from the cluster corrections in Eq. (13) can be trans¬ 
formed into explicit expressions for evaluating contri¬ 
butions by representing each ta factor diagrammati- 
caiiy_, 28 j^ Specifically, the diagrammatic representation 
depicts a bond (which is either a non-sticky or a sticky 
bond, corresponding to the factor Xbond or Xgticky) by 
a solid straight line, an interaction between two nearest 
neighbor united atom groups (i.e., the factor CJXpah ) by 
a solid curved line, and a bending constraint (i.e., the 
factor XYbend) by a dashed curved line. Moreover, non¬ 
stickers and stickers are pictured by crosses and solid 
circles when they participate in weak van der Waals in¬ 
teractions. Equation (13) can then be rendered more 
compact using this diagrammatic representation. 
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where S) and D{B,l,b, S) are the combina¬ 

torial and connectivity factors, respectively, for a given 
diagram with B non-sticky bonds, I interaction lines, b 
bending constraint lines, and S sticky bonds. The con¬ 
tributions to the partition function from the sticky in¬ 
teractions that are retained originate from the zeroth- 
order mean-field term those diagrams with one 
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or more sticky bonds that are consistent with the cur¬ 
rent LCT retention of all diagrams with B + I + S < 4: 
bonds and interaction lines. Dudowicz and Freed^i de¬ 
rived all relevant diagrams with sticky bonds for fully 
flexible linear telechelic polymers. The diagrams includ¬ 
ing the bending constraint lines can then be constructed 
by adding dashed curved lines, connecting pairs of con¬ 
secutive non-sticky bond correlation lines in the relevant 
diagrams. Figures 2-4 provide all diagrams with S > 0 
necessary for including the description of chain semiflexi¬ 
bility within the LCT for linear telechelic polymer melts. 


C. Evaluation of the diagrams with semiflexible 
constraints and sticky bonds 

We briefly describe the evaluation for semiflexible, 
sticky diagrams of I, 6, S) and D{B, I, b, S), the fac¬ 
tors that are necessary to calculate the diagrams. The 
combinatorial factor 7 £)(i?, /, b) has the form, 

^D{B,l,b,S) = SDl{B,l,b,S), (15) 

where sd is the symmetry number and 7(5, Z, &, S) is the 
number of ways of selecting the set of non-sticky bonds, 
sticky bonds, and uncorrelated interacting united atom 
groups from all polymer chains in the system. The sym¬ 
metry number sd for a diagram with S sticky and B 
non-sticky bonds is identical to that for the correspond¬ 
ing diagram with 5-1-5 non-sticky bondsi^ However, 
as illustrated in the example below, certain restrictions 
arise in evaluating 7 ( 5 ,/, 6, 5) for diagrams with 5 > 0 
due to the presence of sticky bondsi^ 

The connectivity factor D{B,l,b, S) for a diagram with 
5 > 0 is identical to that for the corresponding diagram 
where the sticky bonds are replaced by the non-sticky 
bondS)21 and hence, the previous method can directly 
be used to compute D{B,l,b, The connectivity 

factor D{B, I, b, S) can be written as 

5(5, 1, b, S) = (16) 

a 

where a = iV;(iV; — 1) • • • (iVj — Ny + 1) with Ny being 
the number of united atom groups in the given diagram. 
The factor d{B, I, b, S) depends on the lattice structure 
and can be calculated in terms of contracted diagrams?^ 
a set of new diagrams that are obtained by merging sets of 
two or more vertices for the given diagram. Specifically, 
d{B, I, b, S) is evaluated by 

diB,l,b,S)=J2fB,cRB,c. (17) 

C 

where c is a sequential counting index and the co¬ 
efficient fB,c is the product of a contraction factor 

n^i(—l)^^~^(fcA — 1) and the number of ways of form¬ 
ing a contracted diagram with iV' vertices by merging 
the distinguishable vertices fci, ^ 2 ; in the original 

uncontracted diagram. The expressions for the relevant 
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FIG. 5. Illustration of the evaluation of the diagram /) 

( 2 ) 

in Fig. 4. The cumulant diagram /) ' is composed of two 
(2) (2) 

individual diagrams /) { and /) 2 > whose combination yields 
the final result for 


contracted diagrams Rb,c have been provided in Fig. 2 
of Ref. m. 

The diagrammatic representation of the Helmholtz free 
energy fdF = — In IF arises from expanding In IF in a 
Taylor series and from collecting the resulting contribu¬ 
tions into cumulant diagramsi^ Additionally, the Mayer 
/-function is treated using a high temperature expan- 
sioni^ The diagram in Fig. 4 is selected to il¬ 
lustrate the evaluation process for diagrams with sticky 
bonds, and to further clarify the individual contributions 
to 7 d(5, I, b, S) and 5(5, 1, b, S) and the concept of cu¬ 
mulant diagrams. 

('2') 

The cumulant diagram is a combination of two 
individual diagrams, as depicted in Fig. 5. The topol¬ 
ogy of each individual diagram can only be constructed 
in one possible manner by removing the solid curved 
interaction line and labeling all united atom groups as 
distinguishable. Therefore, the symmetry number is 
S£i = 1 for both individual diagrams. Table I pro¬ 
vides a summary for each term of the other factors 
in 7 d( 5,^,6, 5) and 5(5,/, 6, 5). The individual dia¬ 
gram is used as an example to explain the mean¬ 
ing of each term. The bond topology delineated by 
the diagram { means that neither of the two con¬ 
secutive non-sticky bonds can be a bond sequential to 
the sticky bond. This restriction leads to the factor 
H{mN 2 — 2 N 2 e){iTLM — 5){mM — 6)/2\, which is the num¬ 
ber of ways of selecting the combination of one sticky 
bond, two consecutive non-sticky bonds, and uncorre¬ 
lated interacting united atom groups from all polymer 
chains in the system. The symbol A^ 2 e is defined as equal 
to half of the number of runs of two consecutive bonds in 
a single chain for which one of the bonds links a sticker 
with a non-sticker. Hence, N 2 e = 1 for linear chains. The 
diagram i contains seven united atom groups, and thus 
the factor a simply reads a = — i). The factor 

d = —254^2 — 854,5 + 1654,10 + 3254.n with the expres¬ 
sions of 5 b, c given in Fig. 2 of Ref. [H is obtained from 
the contracted diagrams generated when both the inter¬ 
action line and the sticky bond are replaced by non-sticky 
bonds. Combining the results from both individual di¬ 
agrams according to the rules given in Fig. 5 leads to 
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TABLE 1. Expressions for the terms that are necessary for the evaluation of each individual diagram and hence the cumulant 
diagram in Fig. 4. 


Diagram 

a 

d 

7d 

/7,1 

fC2) 

Jr,2 

UloiNi-i) 

nLo(^* “ *) nLo(-^* - *) 

— 2_R4,2 — 8-R4,5 + 16i?4,10 + 32_R4 ,ii 
— 2_R4,2 

H{mN 2 - 2N2e){mM - 5){mM - 6)/2! 
H{mN 2 ){mM - 2){mM - 3)/2! 


the contribution to the free energy from the cumulant 
diagram as 

f(2) ^Kye{^U 24 >^ - - 2N2e4>^)j (18) 

z 

where the variable y = H/Ni arises from the factor 
')D[B,l,b,S). 


D. Contributions to the free energy arising from the 
sticky interactions 


Combining all contributions arising from the sticky in¬ 
teractions [i.e., the zeroth-order mean-field term = 

— 1/Ni InkEJ"^ and the diagrams with S' > 0 in Eq. (20)] 
leads to the final expression for /g, 

4 

/3fs = /3fr^ -Y.Y,y\ (19) 

2=1 


The zeroth-order mean-field term /3/J"-^ reads^l 


=—^a;ln(^a;) -I- (^a; — 2y) \n{(j)x — 2y) 


+y 


l + lnl^ 


- /3es 


( 20 ) 


where x = 2/M is the fraction of stickers in a single chain. 
The coefficients li (i = 1,..., 4) are obtained by collecting 
terms corresponding to a given power of y and can be 
organized in powers of the polymer volume fraction (/), 


j=0 


( 21 ) 


where jmax = 5, 4, 2, and 0 for i = 1, 2, 3, and 4, respec¬ 
tively, and the explicit expressions for Yij are provided 
in Appendix B. 

The variable y in Eqs. (19) and (20) is determined by 
the maximum term method, i.e., by applying the condi- 
tion)21 


dWs) 

dy 


= 0 , 


which, in turn, takes the form. 


{(j)X 




( 22 ) 


0(23) 


TABLE 11. Summary of the meanings of the parameters that 
appear in the free energy / for semiflexible linear telechelic 
melts. 


Symbol 

Meaning 

d or z/2 

spatial dimension 

T 

absolute temperature 

0 

polymer volume fraction 

M 

molecular weight 

Ui (i = 1, ...,4) 

counting indices 

y* 

concentration of the sticky bonds 

e 

van der Waals interaction energy 

ts 

sticky interaction energy 

Ei 

bending energy 


Equation (23) must be solved numerically, producing the 
concentration y* of the sticky bonds that, when substi¬ 
tuted into Eqs. (19) and (20), determines the contri¬ 
butions to the free energy /g arising from the sticky in¬ 
teractions. Notice that y* depends on all molecular and 
thermodynamic parameters including T, (/>, M, e, Eb, and 
Eg. Because the chain’s ends are assumed to be mono¬ 
functional in the present model, the volume fraction of 
the active stickers (i.e., those participating in sticky in¬ 
teractions) is simply 2y*, while the upper limit for y* is 
just y/nax = While the current version of the LCT 

provides no information concerning the concentration of 
sticky bonds in the cyclic clusters, cyclic clusters may 
form, as introduced in Sec. II. 

Finally, the specific free energy / for a melt of semi¬ 
flexible linear telechelic chains appears as 


Pf =13fo - (l)x\n{(j)x) -f {(j)x - 2y*) ln{(j)x - 2y*) 

4 


+ 2 /* 


1 -I- In ( ^ ) -I- /3eg 


-Y.Y,{y*)\ (24) 


where the free energy /o of the reference system is given 
in Subsection III A. For convenience in understanding 
and usage of the theory. Table II summarizes the mean¬ 
ings of the parameters that appear in /. 


IV. DISCUSSION 

The increasing scientific interest in telechelics and 
their technological importance has motivated a num¬ 
ber of theoretical investigations on the self-assembly of 
telechelic polymers. While theories of self-assembly in 
telechelic polymers traditionally employ highly coarse 
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grained models that represent the assembling molecular 
species as a structureless entitythe LCT of Dudowicz 
and Freed^i aims at providing a deep understanding of 
the relation between molecular structure and thermody¬ 
namic properties in telechelic polymers. Therefore, the 
LCT offers the possibility to yield a framework for ex¬ 
tracting and organizing information that is useful for de¬ 
signing telechelic polymer materials. 

The initial study by Dudowicz and Freed^l considers 
for simplicity the fully flexible linear polymers. How¬ 
ever, chain stiffness is known to significantly influence the 
thermodynamics of polymer systems and hence, the role 
of chain semiflexibility in the self-assembly of telechelic 
polymers remains to be investigated within the LCT. 
Moreover, the model of fully flexible linear telechelic 
polymers is not suitable for exploring glass formation 
within the GET— because the characteristic glassy be¬ 
havior appears in the GET only if the polymer chains 
are modeled as being semiflexible. To address the above 
needs, we further extend the LCT for the thermodynam¬ 
ics of linear telechelic polymer melts to include a descrip¬ 
tion of chain semiflexibility. Following the previous treat¬ 
ment of chain semiflexibility within the LCT,— chain 
semiflexibility is treated by introducing a bending energy 
penalty whenever a pair of consecutive bonds from a sin¬ 
gle chain lies along orthogonal directions. We provide an 
analytical expression for the Helmholtz free energy for 
the model of semiflexible linear telechelic polymer melts. 

The present extension provides a theoretical tool for 
investigating the influence of chain stiffness on the 
thermodynamic properties of self-assembling telechelic 
polymers. By combining the present extension of the 
LCT with the AG relation,— we can provide a sim¬ 
ilar generalization of the GET of polymer glass forma¬ 
tion^ thereby permitting exploring the influence of self- 
assembly on glass formation in telechelic polymers. Be¬ 
cause of the great algebraic complexity of the exten¬ 
sion, computations for the influence of chain stiffness on 
the basic thermodynamic properties of self-assembling 
telechelic polymers will be provided in paper H of this 
series.— 
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Cs = + (3ui/2 - U 2 g - 10ui)(/3e)^,(A5) 

Ce = (A6) 

where g = Zpexp(—f]Eb)/zb is called the bending en¬ 
ergy factor, e is the nearest neighbor van der Waals in¬ 
teraction energy, and the counting indices Ui (i = 1,..., 4) 
designate the ratio Ni/M with Ni being the number of 
runs of i sequential bonds in a single chain. The relation 
between Ni and M simplifies to Ni = M — i for linear 
chains. 
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Appendix A: Summary of the coefficients that appear in 
the free energy for the reference system 

The explicit expressions for the coefficients that ap¬ 
pear in the free energy /o for a melt of semiflexible linear 
chains are provided in the following: 


Appendix B: Summary of the coefficients that appear in 
contributions to the free energy arising from the sticky 
interactions 

The explicit expressions for the coefficients that ap¬ 
pear in contributions fs to the free energy arising from 
the sticky interactions for a melt of semiflexible linear 
telechelic chains are provided in the following: 
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where iVse is defined by half of the number of runs of 
three consecutive bonds in a single chain, where one of 
the bonds links a sticker with a non-sticker, and, hence, 
A^ 3 e = 1 for linear chains. 
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